Mon. Not. R. Astron. Soc. 000, [TJED (2012) 



Printed 14 January 2013 



(MN style file v2.2) 



Numerical estimates of the accretion rate onto 
intermediate-mass black holes 



C. Pepe 1 ' 2 * and L. J. Pellizza 1 ' 2 

1 Instituto de Astronomia y Fisica del Espacio, Casilla de Correos 67, Sue. 28, 1428, Buenos Aires, Argentina 

2 Consejo Nacional de Investigaciones CienUficas y Tecnicas, CONICET, Argentina 



Accepted 2013 January 11. Received 2013 January 11; in original form 2012 November 21 



ABSTRACT 

The existence of intermediate-mass (~ 10 3 Af Q ) black holes in the center of globular 
clusters has been suggested by different observations. The X-ray sources observed in 
NGC 6388 and in Gl in M31 could be interpreted as being powered by the accretion 
of matter onto such objects. In this work we explore a scenario in which the black hole 
accretes from the cluster interstellar medium, which is generated by the mass loss of 
the red giants in the cluster. We estimate the accretion rate onto the black hole and 
compare it to the values obtained via the traditional Bondi-Hoyle model. Our results 
show that the accretion rate is no longer solely defined by the black hole mass and 
the ambient parameters but also by the host cluster itself. Furthermore, wc find that 
the more massive globular clusters with large stellar velocity dispersion are the best 
candidates in which accretion onto IMBHs could be detected. 
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1 INTRODUCTION 

The existence of intermediate- mass black holes (10 2 -10 4 M o , 
IMBHs) was suggested, among other observations, by the 
detecti on of ultraluminous X-ray sources in nearby galax- 
ies (see lFeng fc Soriall201ll . for a review). The X-ray fluxes 
and the distances measured for these sources imply lumi- 
nosities (assuming isotropic emission) L > 10 39 erg s" 1 , in 
excess of the Eddington luminosity for an accreting stellar- 
mass compact object [M ~ IQMq). The simplest explana- 
tion for these high luminosities is the presence of accreting 
objects with masses in the range of those of IMBHs, but 
other hypotheses have been proposed. iKording et al.l (|2002l ) 
suggested that the emission of ULXs could be beamed, 
hence impl y ing lo wer luminosities for these sources, while 
iKing et al.l (|200ll ) proposed that ULXs emit indeed in a 
super-Eddington regime with mild geometrical collimation 
of their photon emission (a factor < 10) due to outflows. 
Although in these alternative models accreting stellar-mass 
black holes could explain the observed luminosities, they 
fail to account for sources with L > 10 41 erg s _1 , which 
still require the acc reting compact object to be an IMBH 
jFeng fc Soriall201ll ). 

On the other hand, different theoretical models f or the 
origin of IMBHs have been proposed. According to iFrverl 
(2001), IMBHs would be the dead fossils of primordial 
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(Population III) stars, while iMiller fc Hamilton! ((2002) have 
shown that black holes of 1O 3 M0 could form in globular 
clusters (GCs) as a consecuence of the merger of stellar- 
mass black holes. A s imilar mechanism was proposed by 
iPortegies Zwart et al.l ((2004) , who have shown that a run- 
away collision of massive stars in a GC results in an IMBH at 
its centre. Hence, GCs have become the main candidates to 
host these objects. The extrapolation of the relation between 
the c entral black hole mass and the bulge mass of galaxies 
(e.g.. [Magorrian et al.l(l998l ) also points to the existence of 
IMBHs in GCs. 

Following the predictions of theoretical models, several 
attempts to detect IMBHs at the centres of GCs were done. 
Stellar density profiles and stellar dynamics measurements 
in the central regions of some GCs su ggest the exist ence of 



central objects with masses ~ 10 3 M Q . iMiocchil (|2007l ) devel- 
oped a model for the stellar distribution in GCs, including 
the effects of an IMBH. Those models containing an IMBH 
yield results consistent with the surfa ce brightness and ve - 
locity dispersion profiles obtained by iNovola et al.l (|2006t ) 
for the Galactic GCs NGC 2808, N GC 6388, M80, M13 
M62 a nd M54, and also Gl in M31. IVan den Bosch et all 
(2006) constructed dynamical models of M15 and esti- 
mated the central mass in 3400 Mq. However, the na- 
ture of this concentrated mass can not be distinguished: 
it could be either an IMBH or a large number of stellar- 
mass comp act objects, or a combi nation of both. On the 
other hand, iMcLaughlin et all (|2006l ) fitted the proper mo- 
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tion dispersion profile of 47 Tuc, obtaining an estimate of 
its central point mass o f 100 — 1500 Mq at 68% confi- 
dence level. iNovola et ail l|2010h analyzed the surface bright- 
ness and velocity dispersion profiles of us Cen. Both pro- 
files show a central cusp, a 4000Afp ) IMBH being t h e best 
explanation for these observations. iGebhardt et al.l (|2002t ) 
reported a 2 x lO 4 Af0 point mass at the center of Gl, 
based on photometric and kinematical data in the central 
areas. However, oth er authors arrived at different co nclu- 
sions. For examp le, van der Marel fc Anderson! (|2010h and 
iBaumgardt et all |2003l ') constructed dynamical models for 
oj Cen and Gl, respectively, claiming that the presence of 
an IMBH is not needed to fit the available data. Because of 
this lack of consensus, the study of IMBHs is still an open 
issue and efforts need to be made in order to clarify the sub- 
ject. Moreover, although kinematical and photometric data 
are suitable to show the presence of central mass concen- 
trations in GCs, at present they are unable to determine 
the nature of these concentrations, namely if they comprise 
a single massive object or a collection of stellar-mass rem- 
nants. Complementary data are needed to prove any of these 
hypotheses. 

As in other systems containing compact objects, the de- 
tection of ongoing accretion may help to place constraints on 
the existence and properties of IMBHs. Hence, searches of 
central X-ray sources have been performed in several GCs. 
The detection of central sources with properties consistent 
with those expected from an IMBH has been reported only 
for two clusters: NGC 6388 and Gl (|Nucita et all 120081 ; 
iKong et al.ll2O10l ). In other clusters, only upper limits for 
the X-ray luminosity of a central source have been ob- 
tained, constraini ng it to be lower than ~ ^g3i-32 er g g -i 
l|Maccaronel |2004| . and refereces therein). Assuming that 
the IMBH accretes from the intracluster medium (ICM), 
that the IC M has a density similar to that measured by 
iFreird |200ll ') in NGC 104, and that the radiative efficiency 
is similar t o tha t observed in other black hole systems, 
iMaccaronel |2004l 'l concludes that the accretion rate must 
be far lower than that predicted b y the standard Bondi- 
Hoyle model (|Bondi fc Hovldll944l ) in order to match the 
non-detections of central sources. Based on the correlation 
between X-r ay and radio lumin osity of accreting stellar- mass 
black holes l|Gallo et al.ll2003l ). this author also argues that 
the radio emission of the system would be more easily de- 
tectable than the X-ray emission. However, no central ra- 
dio source has been d etecte d in Galactic GCs which, ac- 
cording to IMaccaronel (|2004h . imposes stron g constraints to 
the m asses of the hypothetical IMBHs. Only Ulvesta d et al.l 
|2007l ) reported the detection of a radio source at about 

1 arcsec from the centre of Gl in M31. They concluded that 
this radio emission is consistent with that expected for a 

2 x 10 4 M© IMBH accreting at the center of the cluster. 
Neverthele ss, this result has been c hallenged by recent ob- 
servations (|Miller- Jones et ai1l2012h . 

It is clear that the detection of (or the failure to detect) 
accretion-powered sources in the centres of GCs is crucial 
for the investigation of the existence and nature of IMBHs. 
In this context, the interpretation of both X-ray and ra- 
dio data have been done in the past using simple models. 
In particular, Bondi-Hoyle accretion has been used, which 
assumes a point mass (the IMBH) accretin g from a homoge - 
neous, static medium. In a previous work ( Pepe et al.ll2012l ) 



we studied the accretion of dark matter onto IMBHs. We 
took into account the gravitational pull of the cluster on 
the accreted matter and found that the accretion rate de- 
pends on the cluster mass, unlike the Bondi-Hoyle accre- 
tion rate that scales as the square of the IMBH mass. The 
Bondi-Hoyle accretion rate is retrieved only for ultrarela- 
tivistic fluids, which are not influenced by the presence of 
the cluster. As the ICM is a non-relativistic fluid, we ex- 
pect that its behaviour departs from that predicted by the 
standard Bondi-Hoyle theory. Moreover, the ICM of a GC is 
not a homogeneous, static medium, because it is fed by the 
red giant stars and cleansed by the passag es of the cluster 
through the Galactic plane (|Robertsl P.9881 . The influence 
of these effects on the accretion rate should be quantified 
to make a proper comparison between models and X-ray or 
radio observations. 

In this work we develop a simple numerical model for 
the accretion flow of the ICM onto an IMBH at the centre of 
a GC, including the cluster gravitational pull and the ICM 
sources. We estimate the accretion rate and its dependence 
on both the IMBH mass, and the GC and ICM properties, 
and compare them with those predicted by the standard 
Bondi-Hoyle model. In Sect.[2]we describe the hydrodynami- 
cal equations of the flow and discuss its main characteristics, 
while in Sect. Owe solve these equations for Milky Way GCs 
with different masses of the IMBHs and temperatures of the 
ICM, stressing the differences between our results and those 
obtained with the Bondi-Hoyle model. Finally, in section 2] 
we discuss the consequences of our results on the search for 
IMBHs, and present our conclusions. 



2 THE MODEL 

2.1 Hydrodynamical equations 

In order to understand the accretion process by an IMBH 
in a GC, we study the dynamics of the ICM in the pres- 
ence of the cluster plus IMBH gravitational potential. We 
aim at computing the accretion rate, to compare it with the 
standard Bondi-Hoyle theory. The ICM is generated by the 
mass loss of the red giants of the cluster. Assuming that the 
distribution of these stars follows the stellar mass density of 
the cluster, and that the average mass loss rate is the same 
for all red giants, the rate at which the density of the ICM 
increases due to the injection of matter by these stars, at 
any position within the cluster can be written as 



ap 



(1) 



where p* is the stellar mass density. The right hand side de- 
scribes the gas injection by the stars at a fractional rate a, 
for which theoretical and observational works obtain estima- 



tions in the range 10 14 -10 11 yr 1 (Fusi Pecci fc Renzinil 



19751; IScott fc R ose 1975; iDupree et al . 1994; Maua s et al 



20061 . and references therein). The cluster plus IMBH grav- 



itational field can be described by the model develo ped by 
iMiocchil (|2007l ). This model is basically a lKind l|l966h model 
with the presence of a black hole at the centre of the clus- 
ter, which modifi es the dynam i cs of the stars in the inner 
region. Following IScott fc Rose! (|l975t ). we assume that the 
ICM is an ideal gas, and that its flow can be considered 
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steady, spherically symmetric and isothermal. Under these 
hypotheses, the flow is governed by continuity and Euler's 
equations. The former is 



1 



-(pr u) 



ap 



(2) 



r 2 dr 

where r is the radial coordinate and u and p the velocity 
and density of the flow, respectively. Euler's equation is 



du 
pu — 
dr 



k B Tdp GM(r)p 



/' 



dr 



aup 



(3) 



where G is the gravitational constant, fee is Boltzmann's 
constant, p is the mean molecular mass of the ICM, and 
M(r) the sum of the stellar mass M*(r) inside radius r and 
the central IMBH mass Mbh- It is assumed here that the 
material is injected with null velocity in the flow. 

These equations can be simplified introducing the vari- 
able q — q/a, where q = pur 2 is the bulk flow, giving 



du 
dr 



2c 



dq * 2 
Tr= pr 



GM (r) 



(« 2 



+ c 2 s )r 2 p* 



(4) 



(5) 



where c s = ksTp -1 is the sound speed. It is important to 
note that, with this definition, we can solve the equations 
independently of a. However, a must be known to calculate 
the density and the accretion rate. 

Equation U can be integrated, giving 



9 = / p*r' 2 dr' + q 
Jo 



M*(r) 
4% 



+ 9o, 



(6) 



where the integration constant go is proportional to the 
accretion rate of the black hole. Although the integration 
should be done from the Schwarzschild radius, this radius is 
negligible with respect to all the scales in our model, hence 
we take it as zero. 

To perform the integration, we define adimensional vari- 



ables £ = rr , tp = 



(and, therefore, ip B — u s a 



), u = q{p Q rl)-\ n*(0 = M*(r)(4jrporg)-\ ^bh = 
AfBH(47rpo''o) 1 j an d = + ^bh, where ro is the 

King radius, po is the cluster central density, and a 2 = 
47rGpo?*o/9 i s the velocity dispersion parameter. With these 
definitions, Eqns. [5] and [6] can be rewritten as 



d$ 



^ 2 - ipi 



dui (ip 2 + tp 2 



e ) 



(7) 



(8) 



These equations describe the flow dynamics in our model. 
They represent an extension of the Bondi-Hoyle expressions 
to the case in which the fluid has sources, and a gravitational 
field other than that of the accretor acts on it. These effects 
are represented by the second and the third terms in brack- 
ets in the right-hand side of Eqn. [8] and the first term in 
the right-hand side of Eqn. [7] The boundary conditions and 



integration method for this set of equations are described in 
the following section. 

2.2 Boundary conditions 

In order to integrate equation [8] boundary conditions must 
be set. Far away from the acretting black hole, the velocity 
of the flow must be positive as there is no matter source 
outside the cluster. On the other hand, on the black hole 
surface, matter can only fall inwards as there is no pressure 
gradient that supports the gravitational pull of the black 
hole and the cluster. As a direct consecuence, there exists a 
stagnation radius £ st at which u = 0. Evaluating Eqn. [7] at 
the stagnation radius gives 



n*(U) = 



-wo, 



(9) 



which means that all the matter ejected by the red giants in- 
side f s t is acretted by the IMBH. This is a direct consequence 
of the hypothesis of stationarity. Thus, the stagnation radius 
separates two distinct regions in space: a regime of accretion 
develops in the inner region, and a wind solution exists in 
the outer one. It is important to highlight that, so far, the 
stagnation radius cannot be uniquely defined by these equa- 
tions. For every stagnation radius there exist a solution to 
the equations. However, they are not all physically accept- 
able, as the Rankine-Hugoniot conditions for the continuity 
of the flow at the stagnation radius require that densities at 

both sides of f a t must be equal. 

As in the Bondi-Hoyle problem (e.g.. lFrank et alj|2002h . 
it can be seen from Eqn. |S]that there exist singular values £ s 
(called sonic radii) for which either the velocity equals the 
sound speed in the medium or the aceleration of the flow 
is null. This is the case when the expression in brackets in 
Eqn. [5] cancels, 



2Vs 



1 duj 



9"(C) 



= 0. 



(10) 



For the same reasons of the classical analysis of the Bondi- 
Hoyle problem, only the transonic curves are consistent with 
the properties of an accretion flow. 

Given that £ s t is not known a priori, the problem was 
solved numerically for a grid of values of £ st between and 
the cluster tidal radius. For each one of these values, sonic 
radii in the inner and outer regions were searched for, using 
a bisection scheme to find the roots of Eqn. 1101 Then, Eqn.[S] 
was integrated (via a fourth order Runge-Kutta scheme) in- 
wards from the outer sonic point and outwards from the 
inner one, up to the stagnation radius. The difference be- 
tween the densities at each side of the stagnation point were 
calculated for each £ st , and interpolated to zero to find the 
true stagnation radius. Once this value was found, we inte- 
grated Eqns. and [8] to obtain the true density and velocity 
profiles. 



3 RESULTS 

3.1 Individual clusters 

We first explore the behaviour of the accetion flow in a given 
GC. The cluster is represented by the value of the concen- 
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Figure 1. Stagnation radius and accretion rate vs gas temperature for different clusters (squares for M 15, triangles for Liller 1, plus 
signs for u> Cen and stars for M 28). The adimensional black hole mass Qbh is ~0.007 in all cases. A decreasing profile arises and the 
curves are placed in the graph ordered with increasing central cluster potential (from bottom to top). 
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Figure 2. Stagnation radius and accretion rate vs black hole mass for different clusters (squares for M 15, triangles for Liller 1, plus 
signs for u> Cen and stars for M 28). The adimensional temperature c^/cr 2 is <~0.34 for all the clusters in the HAR regime (upper set of 
curves). For those in a LAR regime (lower set of curves) the adimensional temperature is around 1, except for Liller 1, for which it is 
~ 3. The curves are placed in the graph ordered with increasing central cluster potential (from bottom to top). It can be seen that the 
dependence of the accretion rate on the black hole mass can be neglected compared to the temperature dependence. 



tration cgc (defined as the ratio of the cluster tidal radius to 
ro), and two of the three parameters ro, po, o. These param- 
ete rs, together wit h the IMBH mass allow us to construct 
the lMiocchil l|2007l) model for the stellar mass M(r) of the 
cluster uniquely. The concentration and IMBH mass define 
the shape of M(r), while the other parameters merely set 
its physical scale. Note that the same is true for our hydro- 
dynamical model, as the flow can be fully expressed using 
non-dimensional variables defined in terms of those scaling 
parameters. Thus, for a given GC there are only two free 
parameters in our model: the non-dimensional temperature 
(or sound speed c s /<t) and IMBH mass Obh. 

With the above considerations in mind, we took four 
sample GCs: NGC 7078 (M15), Liller 1, NGC 6626 (M28) 
and NGC 5139 (uj Cen) , which span the concentration range 
of Milky Way GCs (see Table [1]) , and for each one we con- 



structed different flow models, performing a scan of the two 
free parameters. Although NGC 7078 and Liller 1 have al- 
most the same value for the concentration, they differ in the 
scaling parameters, which allows us to extend the range of 
the non-dimensional free parameters. The temperature cov- 
ers the ra nge 5000—15000 K, which is the range expected for 
the I CM i|Scott fe Roselll975l ; lKnapplll996l ; [Priestley et ail 
12011]) . The IMBH mass ranges from 10 2 M© up to 1O 4 M . 
Heavier masses would produce a strong effect on the stellar 
dynamics and structure of the cluster, which is not observed, 
while for lower masses the assumption of the IMBH at rest 
at the center of the cluster would not hold. 

In Figs. [1] and [2] we show the stagnation radius £ s t 
and the (non-dimensional) accretion rate f2(£ s t) as a func- 
tion of the free parameters c s /<t and Qbh, respectively, 
with the other parameter kept fixed. It is worth pointing 
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Cluster 


<7 


ro 


P0 




CGC 




(km s — 1 ) 


(PC) 


(M Q 


pc" 3 ) 




NGC 7078 (M 15) 


10.8 


0.42 


1.12 


X 


10 5 


2.29 


Lillcr 1 


5.05 


0.15 


1.90 


X 


10 5 


2.3 


NGC 5139 (w Cen) 


10.2 


3.58 


1.41 


X 


10 3 


1.31 


NGC 6626 (M 28) 


7.8 


0.38 


7.24 


X 


10 4 


1.67 



Table 1. Globular cluster parameters taken from [Harri 

The central density was estimated assuming a mass-luminosity 

relation Mq ~ Lq. 
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out that the true accretion rate M can be recovered via 
M = af2(£ s t)po'"o, with a suitable value for a (see Sect. [4] 
for a discussion) . It can be seen from the left panel of Fig. [1] 
that the stagnation radius decreases with the ratio c s /a. 
This can be easily understood in terms of the energetics of 
the gas: for higher temperatures the gas has more energy and 
can escape from inner regions of the cluster, where the grav- 
itational well is deeper, hence moving the stagnation radius 
inwards. This results in a lower accretion rate (right panel). 
It is interesting to note that the stagnation radius curve 
steepens at r st ~ ro, where c s /a ~ 1. The rapid increase of 
the stellar mass of the cluster at this radius results then in 
an abrupt change in the accretion rate. This abrupt change 
suggests the existence of two accretion regimes: at high tem- 
peratures (cs/a > 1), the stagnation radius is located in the 
central region of the cluster (r s t < ro) resulting in a low 
accretion rate (LAR), while at low temperatures (c s /a < 1) 
the stagnation radius is far from the centre (r st 3> ro) and 
the accretion rate is high (HAR). Note that the accretion 
rate differs by almost three orders of magnitude between 
the LAR and the HAR. Fig. [T] shows another characteristic 
differentiating the HAR from the LAR. In the former, the 
stagnation radius and the accretion rate increase with the 
cluster concentration cgc- 

Fig. [2] shows the dependence of the stagnation radius 
and accretion rate on the non-dimensional IMBH mass, at 
fixed non-dimensional temperature. Two temperatures were 
explored, one resulting in the LAR and the other in the 
HAR. Both the stagnation radius and the accretion rate are 
almost independent of Q,bh in the HAR regime. This hap- 
pens because far from the centre the IMBH has no influence 
on the stellar distribution, and the stellar mass increases 
very slowly. On the other hand, in the LAR regime there is 
a strong dependence on the IMBH mass, because the stag- 
nation radius is near the sphere of influence where the stellar 
distribution becomes dominated by the IMBH. Note that in 
this regime, the accretion rate varies also with the cluster 
concentration, decreasing as the latter increases. In Fig. [3] 
we show the ratio between the stagnation radius and the 
accretion radius defined as r acc = GMbh/c;?. It can be seen 
that in the HAR regime, since the stagnation radius remains 
independent of the black hole mass, the ratio r st /r acc de- 
creases with increasing Mbh- In the LAR regime, this ratio 
remains almost constant implying that M scales as ~ Mbh 
like in the traditional Bondi-Hoyle model. However, the ra- 
tio is 1-2 orders of magnitude greater than unity, and hence 
the accretion rate at a fixed IMBH mass is much higher in 
our models than in the Bondi-Hoyle case. 



100 1000 10000 

Figure 3. Ratio between the stagnation radius and the accretion 
radius of black hole. In the LAR regime the stagnation radius 
behaves like the accretion radius. 



3.2 Milky Way globular clusters 

In the previous section, we have shown that in some cases 
there is a dependence of the accretion rate on the cluster 
properties. Here we explore this dependence by applying 
the scheme detailed in Sect. [2] to the Milky Way G Cs with 
well d etermined parameters (listed in the catalogue of lHarrisl 
Il99rj) . to construct flow models for different temperatures 
and IMBH masses. It is worth pointing out that the results 
from Sect. [5] extend to all the clusters in the catalogue and 
the same reasoning can be applied to them. For each model, 
the accretion rate was computed using the same value of a 
as in the previous section. In some cases, the models were 
discarded because the IMBH mass was not much lower than 
the cluster mass, while in others no stagnation radius was 
found for the hottest temperatures (the ICM escapes com- 
pletely as a wind). We searched for correlations between the 
accretion rate and the properties of the cluster, defined by 
their scaling parameters and masses. We found no trend for 
the accretion rate M with ro or po, but a clear correlation 
with the cluster mass Mgc and velocity dispersion parame- 
ter a. 

We show in Fig. [5] the dependence between M and a 
for a conservative value of the IMBH mass, Mbh = lOOOM© 
and three different temperatures T = 5000, 10000, 12600 K 
(left panel), and for a fixed gas temperature (T = 5000K) 
and different IMBH masses M BH = 1000, 4000, 10000M Q 
(right panel). The clear trend observed in the left panel is 
understood in terms of the results of the previous section. 
For a fixed value of T, an increase in a implies a decrease 
in Cs/a, changing the flow from the LAR upwards to the 
HAR in the curves defined in Figure [T] This explains why 
the GCs with higher a are more likely to be found in a HAR 
regime. The change in the location of the curve as the tem- 
perature increases is also consistent with this interpretation. 
The dispersion of this correlation is in part due to the differ- 
ent concentrations of the GCs, and in part to the fact that 
M scales as poro, which is different for each cluster. The 
right panel of Fig. [4] shows that the effect of the IMBH mass 
for a fixed temperature is small. 

A trend of increasing accretion rate with the cluster 
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mass can also be seen in Fig. [5] for high accretion rates. 
This is explained by the fact that these clusters are accret- 
ing in the HAR regime, for which the stagnation radius is 
well outside the cluster core and hence encloses almost the 
whole cluster stellar mass. As M is proportional to the stel- 
lar mass inside r st , the correlation with Mgc arises. In other 
words, most of the mass ejected by the stars in the cluster 
is accreted by the IMBH. As the cluster is more massive, 
more mass is available to be accreted, which explains the 
clear trend with Mgc in the HAR regime. Note that this 
trend is independent of the temperature and IMBH mass 
(as far as the accretion proceeds in the HAR regime), and it 
is very tight because both M and A^gc scale with the same 
combination of parameters of the cluster (po^o), reducing 
the dispersion due to the scaling from non-dimensional to 
physical variables. For clusters in the LAR regime there is 
no clear trend, as the small stagnation radius decouples the 
accretion rate from the cluster mass. According to the re- 
sults of this section, the clusters with large values of a and 
high Mgc are most likely to be in a HAR regime, and are 
therefore the best candidates to perform detections of the 
accretion onto an IMBH. 



3.3 Differences with the Bondi-Hoyle model 

The model presented in this paper differs from that of 
iBondi fc Hovlel | 1944 ) in some very fundamental aspects. 
Our model takes into account the gravitational potential 
of the cluster, which differs from cluster to cluster, and 
also includes the constant injection of gas into the ICM 
by the red giants of the cluster, hence the black hole is no 
longer accreting from a static and infinite medium. These 
two features of our model produce a very important conse- 
quence, which is its main difference with the Bondi-Hoyle 
result: the accretion rate not only depends on the black 
hole mass and gas temperature (or sound speed), but also 
on the cluster properties. Figure [S] shows the comparison 
between our model (for different clusters; all of them at 
T = 10000 K), the Bondi-Hoyle accretion rate, and the 
Eddington rate Msdd = -^Edd/c 2 , where c is the speed of 
light and L Ed d = 1.26 x 10 38 (M BH /M Q ) erg s 1 is the Ed- 
dington luminosity of the IMBH. The Bondi-Hoyle accretion 
rate was calculated as Mbh = 4-kG 2 M% h p a c~ s , where p a is 
the ambient gas density. We used the value p a = 0.2 cm -3 
(|Freire|[200lh . which is the same value used by other authors 
(e.g., iMaccarond l2004h to compute IMBH X-ray luminosi- 
ties. To make the results of our models comparable to those 
obtained with the Bondi-Hoyle scenario, for each of our mod- 
els we have chosen the fractional mass loss rate of the stars 
a so that the density at the stagnation radius matches p a . 
This choice is justified because it makes both models to have 
the same density at the point where the fluid is at rest. The 
values of a obtained are in the range 10 -14 — 10~ n yr _1 , in 
rough agreement with (but somewhat lower than) the few 
observational estimations of this parameter (|Scott fc Rose! 
ll975l ; |Priestlev et alJlioTlh . 

It can be seen from Fig. that the accretion rate in 
the HAR regime (such as that of NGC 7078 in this plot) 
no longer scales as M| H , and that it is cluster-dependent. 
However, in the LAR regime the behaviour is similar to that 
of Bondi-Hoyle models, although the absolute value of M 
is one order of magnitude higher. Another interesting re- 
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Figure 6. Comparison between our accretion rate and Bondi- 
Hoyle accretion rate. We also show the Eddington accretion rate 
for a reference. 



suit seen in this figure in that in the HAR regime, IMBHs 
reach accretion rates as high as the Eddington values (as- 
suming a 100% efficient conversion of gravitational energy 
into radiation). This result is due to the development of a 
cluster- wide flow in the HAR regime, which carries matter 
from the outer regions onto the IMBH to produce a high 
accretion rate. It cannot be reproduced by the Bondi-Hoyle 
model unless a very high density or very low temperature 
are assumed. This result is important for the explanation 
of the emission of ULXs, as will be discused in the next 
section. It is worth pointing out that these flows with super- 
Eddington accretion rates do not imply super- Eddington lu- 
minosities as the radiation efficiency is usually much lower 
that 1. Hence, these flows can still proceed without being 
stopped by radiation pressure. 



4 DISCUSSION 

In this work, we developed a simple model for the accre- 
tion of the ICM of a globular cluster onto an IMBH at its 
center. Bearing in mind that the detection of ongoing accre- 
tion would be a strong piece of evidence for the existence 
of these elusive black holes, our aim was to refine the pre- 
dictions of the expected accretion rate, improving on the 
(usually assumed) classical Bondi-Hoyle accretion. In par- 
ticular, we explored the consequences of including the effect 
of the gravitational field of the cluster and the fluid sources 
on the flow. Our models aims at assessing the influence of 
the cluster as a whole on the accretion flow, which was sug- 
gested to be impor tant in the case of cosmological fluids such 
as dark matter bv lPepe et all (|2012T ). 

We made several assumptions in our model, namely that 
the flow is stationary, spherically symmetric, and isothermal. 
While the first two can be regarded as working hypotheses 
made to simplify the computations, th e latter can be a goo d 
approximation for the whole cluster (|Scott fc Rose! Il975l ). 
However, it may break down very near the black hole, where 
the transformation of gravitational into thermal energy is 
higher and may heat the gas in timescales shorter than the 
cooling one. Although our model could not describe the dy- 
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clusters with low <r. The clusters located in the right end of the curve are those clusters accreting at a HAR regime. The dotted lines 
indicate the a values corresponding to the clusters from Sect. [2] and the particular case of NGC 6388 discussed in Sect. l4l 




A 

□ 



u Ie-07 

D 

O 

< 

le-0S 




+4*1+ 



A^ 



Figure 5. Accretion rate vs globular cluster mass for three different temperatures (left panel): 5000K (triangles), 10000K (plus signs) 
and 12300 K (squares). Accretion rate vs globular cluster mass for three diffcrente black hole masses (right panel): IOOOMq (triangles) 
, 4OOOM0 (plus signs) and 10000Mq (squares). A soft trend can be observed, although there is a wide spread of the data. The dotted 
lines indicate the Mqq values corresponding to the clusters from Sect. [2] and the particular case of NGC 6388 discussed in Sect. [4] 



namics of the flow in these spatial scales, it is still useful to 
assess the existence of cluster-wide flows, whose dynamics 
is governed by phenomena occurring at much larger spatial 
scales, and estimates the order of magnitude of the accre- 
tion rate produced by these flows. To analyse the details of 
the flow very near the black hole, we are currently develop- 
ing more complex models which include several heating and 
cooling mechanisms. 

Our main results are the following. First, large scale 
flows can form in globular clusters with IMBHs, due to the 
influence of the gravitational field of the cluster and the 
mass loss of the cluster stars. These flows show a stagnation 
radius at which the velocity is null. Outside this radius, the 
ICM escapes from the cluster as a wind, as the energy of the 
flow allows it to overcome the cluster potential well. Inside 
the stagnation radius the ICM is retained by the potential 
well, and an accretion flow onto the IMBH develops. As our 



models are stationary, the accretion rate onto the IMBH is 
determined by the stagnation radius and the mass-loss rate 
of the stars. 

The accretion rate onto the IMBH in a given cluster 
depends mainly on the ratio of the sound speed in the ICM 
to the stellar velocity dispersion parameter. This is a con- 
sequence of the energet i cs of the flow, and the same effect 
observed bv lPepe et al.l |2012l ) for dark matter. The sound 
speed (related to the temperature of the flow) measures the 
internal energy of the ICM, while the stellar velocity disper- 
sion parameter indicates the strength of the cluster gravita- 
tional field. As the ratio of these two magnitudes increases, 
the ratio of the internal to gravitational energy of the flow 
increases as well, making easier for the ICM to escape from 
the cluster. Therefore the stagnation radius decreases, the 
wind becomes stronger and the accretion flow weaker, di- 
minishing the accretion rate. 
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The stellar mass distribution of globular clusters has 
also important effects on the flow. These clusters show a 
large concentration of their mass within a few core radii of 
their centres, while the potential well of the cluster extends 
to far beyond (typically to tidal radii one or two order of 
magnitudes larger than the core radius) . This translates into 
a steep increase of the accretion rate when the temperature 
is low enough for the stagnation radius to reach the core 
radius, as most of the cluster stellar mass (and hence most 
of the stellar mass loss) is enclosed within the latter. This 
steep increase separates two accretion regimes with different 
properties, one with a high accretion rate at low tempera- 
tures (HAR), and the other with a low accretion rate at 
higher temperatures (LAR). 

In the HAR regime, for a fixed temperature the ac- 
cretion rate is high, proportional to the cluster mass, and 
almost independent of the IMBH mass, as far as the IMBH 
is not massive enough to severely change the whole cluster 
mass distribution (however, if this were the case, strong sig- 
natures of the presence of the black hole should be found in 
the stellar distribution and dynamics). The dependence of 
the accretion rate on the cluster mass instead on the IMBH 
mass is due to the fact that in the HAR regime, the stagna- 
tion radius is in the outer region of the GC, and the IMBH 
accretes a major fraction of the mass lost by the stars, which 
is proportional to the cluster mass. Note the difference with 
the classical Bondi-Hoyle model, in which the accretion rate 
scales as the square of the black hole mass. 

The LAR regime is qualitatively different, showing ac- 
cretion rates several orders of magnitude lower, which de- 
pend strongly on the IMBH mass. The accretion rates in 
the LAR regime are still one order of magnitude greater 
that the Bondi-Hoyle prediction for similar boundary condi- 
tions, and depend on the IMBH mass because the stagnation 
radius is near or inside the IMBH sphere of influence. This 
dependence varies from cluster to cluster due to the differ- 
ent stellar-mass profiles of the clusters, but on average it 
mimicks the Mg H dependence of the Bondi-Hoyle model. 

Another strong assumption of our models is that the 
IMBH is at rest at the centre of the cluster. This approxi- 
mation holds because the mass of the IMBH is far greater 
that the mass of the stars in the cluster. Assuming that the 
stars in the cluster have a mean mass of ~ 0.5Mq and a 
mean velocity of the order of a, energy equipartition im- 
plies that the IMBH would have a velocity of the order of 
a /2^J Mbr/Mq. At these velocities, and taking into account 
that c s /a ~ 1 in our models, the typical correction to the 
accretion r ate due to the motion of the IMBH with respect 
to the flow (|Hovle fc Lvttletonll 19391 ) would be less than 1%. 
Hence, neglecting this effect is justified for these models. 

The application of our model to the Milky Way globular 
clusters has shown that the higher accretion rates are pre- 
dicted for those cluster with the largest values of the velocity 
dispersion parameter a. However, the higher the tempera- 
ture of the gas, the higher the minimum value of a for GCs 
accreting at high rates. This result suggests that the sig- 
natures of the accretion onto IMBHs must be searched for 
in globular clusters with high velocity dispersion parameters 
and low gas temperatures. All the clusters that were observa- 
tionally proven so far for the existence of IMBHs satisfy the 
first criterion. The satisfaction of the criterion on gas tem- 
perature is difficult to assess, as measurements of the ICM 



thcrmodynami cal state are s t ill rar e (e.g., lFreirdl200ll ). As 
pointed out bv lScott fc Rose] (|l975l ), the temperature of the 
flow depends mainly on that of the radiation field to which it 
is exposed, hence a proxy for the former would be the num- 
ber of UV sources within the cluster, such as blue horizontal- 
branch stars or blue stragglers. These sources would heat 
the gas, preventing the accretion flow to develop out to large 
scales, and hen ce decreasing the accretion rate. Interestingly, 
iMiocchil (|2007T ) argues that extreme horizontal-branch stars 
with strong UV fluxes might be the result of the stripping 
of normal stars passing near the IMBH. If this is indeed the 
case, the presence of the IMBH would help in heating the 
flow, decreasing the accretion rate. 

From the results of our model we can estimate the X-ray 
luminosity Lx due to the accretion process onto the IMBH. 
However, in our model the luminosity depends linearly on 
two poorly constrained parameters, the fractional mass loss 
rate of the stars a and the accretion efficiency of the flow, 
e = Lx/Mc 2 . If we adopt the standard values, e = 0.1 for 
accreting stellar-mass black holes and a — 10 -11 yr _1 , the 
luminosities produced by the accretion flows of our models 
would be in the range 10 37 — 10 41 erg s _1 . The high end 
of this range is in good agreement with the luminosities of 
ULXs. If, as claimed in a few cases, ther e are ULXs position- 
ally coincident with ex tragalactic GCs (|Angelini et al 1 l200ll : 
M accarone et al.ll201ll . and references therein), our models 
suggest that these systems may contain IMBHs accreting 
in the HAR regime at the centres of GCs, with standard 
accretion efficiencies. 

In the case of NGC 6388 our models predict (based on 
the estimation of the accretion rate and assuming the ef- 
ficiencies discussed above) X-ray luminosities in the range 
10 38-40 erg s , at least 5 orders of magnitude higher than 
the observed value of Lx, NGC6388 = 2.7 x 10 33 erg s _1 . 
However, the prediction can be re conciled with obser- 
vations if either a or e are lower. iNucita et all (|2008T l 
reaches the same conclusion about the efficiency. It is 
worth mentioning that different theoretical models have 
been developed for accretion flows with very low efficien- 
cies, such as advection-dominated accretion flows (ADAFs, 
' — — 11 1 

Narayan & Yi 1994) , jet-domi nated accretion flows (JDAFs, 
Fender. Gallo fc Jonkerll2003l), and low-r adiative efficiency 
accretion flows (LRAFs, Qua taert 2001). For example, in 
the case of ADAFs the efficiency scales as e = M /M^dd for 
M < O.lMEdd, giving in our case estimations of the luminos- 
ity 2-3 orders of magnitude lower, approaching the obser- 
vational limits. On the other hand, the value of a is highly 
uncertain, relying in t heoretical models and with few ob- 
servational constraints dScott fc Rose! 19751 ; iMcbonald et al.l 
l201ll ; iPriestlev et~aLl l201lh . Given that the low luminos- 
ity limits of our predicted range correspond to low IMBH 
masses, we conclude that such an IMBH accreting in the 
LAR is a possible explanation for the central engine of the 
central X-ray source in NGC 6388. The case of Gl is slightly 
different, as this cluster is so massive that the luminosity in 
the HAR is Lx ~ 10 42 erg s , and the development of a 
LAR requires temperatures of several times 10000 K. As the 
observed luminosity is only Lx — 2 x 10 36 erg s -1 , it is diffi- 
cult to reconcile it with the predictions, unless the Gl ICM 
is strongly heated, the accretion efficiency is extremely low, 
or its stars lose mass at very low rates. The Galactic GCs 
with undetected central X-ray sources, with measured upper 



© 2012 RAS, MNRAS 000, \T}M 



Accretion rate onto IMBHs 9 



limits of io 31-32 erg s _1 for their X-ray luminosities are in 
the same case. 

The models presented in this work explore the conse- 
quences of taking into account the gravitational field of the 
GC and the mass loss by stars in the computation of the ac- 
cretion rate onto an IMBH in the centre of a GC. Although 
they are very simple, they allow us to get some insight in 
the development of cluster-wide flows that may feed these 
compact objects. Many improvements can be made on the 
models, such as relaxing the isothermal hypothesis to include 
the detailed physics of gas heating and cooling, or including 
the radiation pressure on the flow due to the X-ray emission 
of the flow itself. At present we are working in these tasks, 
which require a more complex numerical approach to the 
problem. We expect that an improved physical model that 
gives more precise predictions on the accretion rates and lu- 
minosities, together with detailed stellar evolution models 
that predict stellar mass loss rates, deeper X-ray measure- 
ments and ICM observations, we would be able to determine 
the existence or not of these elusive compact objects and put 
limits on their masses. 
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